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ABSTRACT 

We report that electron-capture supernovae (ECSNe), arising from collapsing oxygen- neon- 
magnesium cores, are a possible source of 48 Ca, whose origin has remained a long-standing puzzle. 
Our two-dimensional, self-consistent explosion model of an ECSN predicts ejection of neutron-rich 
matter with electron fractions Y c w 0.40-0.42 and relatively low entropies, s « 13-15 fcs per nucleon 
(fee is the Boltzmann constant). Post-processing nucleosynthesis calculations result in appreciable 
production of 48 Ca in such neutron- rich and low-entropy matter during the quasi- nuclear equilibrium 
and subsequent freezeout phases. The amount of ejected Ca can account for that in the solar in- 
ventory when we consider possible uncertainties in the entropies. ECSNe could thus be a site of 48 Ca 
production in addition to a hypothetical, rare class of high-density Type la supernovae. 
Subject headings: nuclear reactions, nucleosynthesis, abundances — stars: abundances — supernovae: 
general 
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1. INTRODUCTION 

The origin of the most neutron-rich (n-rich) stable iso- 
tope among the iron-group nuclei and lighter elements, 
48 C a, remains a long-standing mystery of nuc leosynthesis 
(see lHartmann et al.|[l985l: iMever et al.lll996l for the his- 
torical background). 48 C a constitutes only 0.1 87% of the 
solar calcium abundance (jAsplund et al.l l2009f), but is 47 
times more abundant than its next heavier stable isotope, 
46 Ca. From the point of view of nuclear stability, this 
can be attributed to the doubly magic structure of 48 Ca 
with Z = 20 and N — 28. Thermal equilibrium, such 
as nuclear statistical equilibrium (NSE) or quasi-nuclear 
equilibrium (QSE), thus seems to be responsible for the 
creation of 48 Ca, in which each isotopic abundance is a 
strong function of its binding energy per nucleon. Such a 
formation of 48 Ca demands a highly n-rich environment 
with Y c (number of protons per nucleon) being close to 
t hat characterizing the st ructure of 48 Ca, i^.nuc = 0.417. 

lHartmann et al.l (|1985l) suggested n-rich NSE as the 
production mechanism of 48 Ca, based on the assump- 
tion that such conditions would be met in t he innermost 
eject a of core-collapse supernovae (CCSNe). IMever et al.l 
(1996) argued, however, that QSE was responsible for 
the creation of 48 Ca. They showed that a low-entropy 
environment, characterized by 4> < 1) led to the ap- 
preciable creation of 48 Ca, where (f> = 0.34Xg/p5 is the 
photon-to-nucleon ratio with Tg being the temperature 
in 10 9 K an d ps the mat t er den sity in 10 5 g cm' 3 . For 
this reason, IMever et al.l (|1996D concluded that CCSNe 
were excluded as the origin of 48 Ca because of their 
high entropies (</> > 1) in the neutrino-driven ejecta. 
Such a high 4> leads to an q-rich freezeout from NSE 
(IWooslev fe HoffmanllT99l IMever et al.N1998al) followed 
by QSE, in which nuclear stability prefers heavier nuclei 
than 48 Ca. Instead, they concluded Type la supernovae 
(SNe la) as the unique site that could achieve such low 
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<j) environments. IWooslev! (|1997|) explored SN la models 
with a central density of (4-8) x 10 9 g cm -3 (needed to 
obtain Y c ~ 0.40-0.42), which was substantially higher 
than the typical values of about 2 x 10 9 g cm -3 . He found 
the production of 48 Ca/ 56 Fe to be about 100 times larger 
than its solar ratio. Given the observational fact that 
about one half of th e Galactic iron comes from SNe la 
(|Timmes et al.lll995D . IWooslevl (|1997l) estimated an event 
rate of such high-density explosions no more than a few 
percent of the observed SNe la rate. 

In this Letter, we propose that electron-capture SNe 
(ECSNe), a sub-class of CCSNe arising from collapsing 
O-Ne-Mg cores, can be additional sites of 48 Ca produc- 
tion. A recent two-dimensional (2D), self-consistent sim- 
ulation of an ECSN predicts the ejection of n-rich matter 
of Y c ~ 0.40-0.4 2 with s ~ 14fc B / r iuc (fc B is the Boltz- 
mann constant; Uanka et alj 120081 120121: iWanaio et all 
120111) . These entropies, or <f> ~ 1.4, are only slightly 
above t he limit of the criter ion for 48 Ca formation consid- 
ered bv lMever et~aT1 (|1996l ). We utilize this ECSN model 
(§2) for nucleosynthesis calculations and discuss the pro- 
duction mechanism of 48 Ca (§ 3). The mass-integrated 
nucleosynthetic yields are then compared with the solar 
abundances to test if our ECSN model accounts for the 
48 Ca abundance in the solar inventory (§ 4). 

2. ECSN MODEL 

The nucleosynthesis analysis makes use of about 2000 
representative tracer particles, by which the thermody- 
namic histories of ejecta chunks were followed in our 
self-consistent 2D hydrodynamic calculation of an ECSN 
(| Janka et alj200a,l20TirWariaio et alTMTl . The model 
was computed with a ray-by-ray-plus treatment of the 
energy-dependent neutrino transport ( iRampp fc Jankal 
l200Ct iBuras et al.ll200& IKitaura et all 120061 ). The pre- 
collapse model of the O-Ne-M g core emerged from the 
evolution of an 8.8M Q star (|Nomotol 119871) . In our 
model, a neutrino-powered explosion sets in at about 
100 ms after core b ounce with an energy of about 10 50 erg 
(| Janka et all 12008) . Because of the very steep density 
gradient near the core surface, an interesting amount of 
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Fig. 1. — Top: Ejecta masses vs. Y c at Tg = 5 for our ECSN 
model. The width of the Y e -bins is chosen to be AY C = 0.005. The 
total ejecta mass from the core is 1.14 X 10~ 2 Mq. Also indicated 
by a dashed line is Y e ,nuc = 0.417. Bottom: s and cf> as functions 
of Y e at Tg = 5 for all tracer particles. 

n-rich matter with relatively low entropy (s w 13-15 fee 
for Y e < 0.42) gets ejected (Figure ED topjj This hap- 
pens only in the multi-dimensional case, which allows for 
the rapid expansion of neutrino- heated gas in buoyant, 
mushroom-like plumes, whose expansion is too rapid for 
v e and v e absorptions to lift Y c to values closer to 0.5, 
in contrast to the ID case or more massive Fe-core pro- 
genitors, where the ejecta expansion is much slower and 
neutrino absorption proceeds longer. The entropies are 
lower for the lower Y c particles as a result of less neu- 
trino energy deposition. The photon-to-nucleon ratios 
(Fig. [TJ bottom) are approximately one-tenth of the en- 
tropies, showing however a slightly steeper gradient for 
Y c < 0.42. The expansion timescales of the particles, 
defined as the e-folding time of the temperature drop 
below 0.5 MeV, are r oxp = 50-100 ms and are largely 
independent of Y e . 

3. 48 Ca PRODUCTION 

The nucleosynthetic yields are obtained with t he re- 
actio n network code described in I Wanajo et al.l (|2009l 
1 2 llh by e mploying the latest r eaction library of REA- 
CLIB V2.0 (jCvburt et al.l 120 10ft . Using thermodynamic 
trajectories directly from the ECSN model, the calcu- 
lations are started when the temperature decreases to 



Fig. 2. — Top: Final mass fractions of 48 Ca for the tracer particles 
in the range of Y c < 0.45. Also indicated by a dashed line is 
Kj.nuc = 0.417. The result of the unchanged model is shown in 
red, and those with the densities multiplied by factors of 1.3, 2.0, 
10 and divided by factors of 1.3, 2.0 are given in different colors. 
Bottom: Mass numbers of nuclei that have the greatest production 
factors (X/Xq) as functions of Y e (slightly shifted in the vertical 
direction for visibility). 

Tg = 10, assuming initially free protons and neutrons 
with mass fractions Y^i and 1 — y Cj i, respectively (Y e ^ is 
the value at T 9 = 10). ' 

The final mass fractions of 48 Ca for the tracer parti- 
cles in the range of Y c < 0.45 are shown as a function 
of Y c by red dots in the upper-panel of Figure [2j We 
find appreciable production of 48 Ca in the range below 
its F Cinuc (dashed line), peaking at Y c w 0.4(0. 48 Ca is 
always made as itself, its double-magic nature affording 
a special degree of stability in nuclear equilibrium. 

To examine the nucleosynthesis of 48 Ca in detail, let 
us pick a representative tracer particle with Y c = 0.417, 
s = 14.3/cb/iiuc ((/> = 1.45), and r cxp = 56.7 ms. Note 
that the Y e is slightly greater than Y c \ = 0.414 as a result 
of neutrino reactions and the a-effect ((Mc Laughli n et aTl 
Il996t iMever et al.lH998bl ). The abundances (defined by 
Y = X/A] number per nucleon) of n, p, a, heavy nu- 
clei (Z > 2), and 48 Ca are displayed as a function of 
descending temperature in panel of Figure [3] (upper- 
left). We find that the abundance of heavy nuclei is 
fixed around Tg = 6.24. This is a consequence of the 
fact that the 3-body reactions a(aa,7) 12 C (hereafter 
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result in the drop-off of X( 48 Ca) around Y e < 0.390. 
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Fig. 3. — Abundances of n, p, a, heavy nuclei (Z > 2), and 48 Ca as functions of descending temperature for the tracer particle with 
Y c = 0.417, s = 14.3fc B /nuc (0 = 1.45), and fcxp — 56.7 ms. The result of the unchanged model is shown in the upper-left, and those 
for the densities divided by a factor of 1.3 and multiplied by factors of 2.0, 10 are shown in upper-right, lower-left, and lower-right panels, 
respectively. The long-dashed lines mark the NSE-freezeout temperatures. Also indicated by dashed and long-short-dashed lines are the 
temperatures at which the a-consumption timescales evaluated by 3a and aan, respectively, become longer than r cxp . 



3a) and a(cm,7) 9 Be followed by 9 Be(a, n)C 12 (here- 
after aan) become slow as the matter expands. The 
a-consumption timescale, r^ a = —Y a /{dY a /dt) evalu- 
ated by 3a, first exceeds r exp at Tg = 7.22 (dashed line). 
The gateway from light to heavy nuclei is closed when 
the timescale r aan evaluated by aan as well exceeds r eX p 
(T 9 = 5.41; long-short-dashed line). In general, NSE 
freezes between the 3a and aan freezeouts. We define 
the NSE-freezeout temperature (Tg = 6.24; long-dashed 
line) at which the timescale of heavy abundance forma- 
tion, Thoavy = Yhcuvy / (dYhcavy / dt) , exceeds T CX p. Once 
the NSE freezes, the number of heavy nuclei in QSE 
cannot decreas e because of their n et photodisintegration 
being too slow ()Mever et al.lll996l) . The fixed number of 
heavy nuclei is in fact the key condition of QSE. 

We find that the transition from NSE to QSE takes 
place with the number of a's slightly greater than that 
of heavy nuclei. The abundances become redistributed 
by absorptions of light nuclei decreasing the abundance 
of 48 Ca and increasing the abundances of heavier nuclei 

5 In QSE, individual reactions, which are in general much faster 
than the expansion timescales of matter, are irrelevant. What de- 
termine the abundance distribution are nuclear binding energies 
per nucleon. Individual reacti ons play roles only after a freezeout 
from QSE (Tg m 4; see § 4 in lMever et~allll998al) . Note also that 
a's are in nuclear equilibrium with free nucleons throughout the 
QSE phase. All n, p, a thus play roles for readjustment of the 
QSE abundances. 



in the QSE cluster (Figure HJ. In the late QSE phase 
(Tg w 5.5-4), however, the 48 Ca abundance recovers in 
response to a drop-off of the a abundance. The num- 
ber of charged light nuclei becomes insufficient to trans- 
mute 48 Ca into heavier nuclei. The abundance redistri- 
bution is predominantly due to photodisintegration and 
n-absorption, which cannot make a net upward shift in Z. 
Instead, 48 Ca increases in the local equilibrium around 
its double-magic intersection, lasting until the freezeout 
from QSE at Tg rs 4. In the end, however, the most 
abundant isotope is 66 Zn, the daughter of p-magic 66 Ni 

(Fig. ED. 

How can we keep 48 Ca increasing in the QSE phase? 
Obviously, the presence of a's is the cause of shifting light 
QSE abundances toward a greater Z in the above case. 
We thus test the sensitivity of the 48 Ca production to 
entropy, the key quantity that controls the a abundance. 
Figure [3] gives the results in which densities are divided 
by 1.3 (upper-right) or multiplied by 2.0 (lower-left), 10 
(lower-right) throughout the calculations for the same 
tracer particle. The 0's (and roughly s's) decrease or in- 
crease by the inverse of the density scaling factor in each 
case. We find that the "p/1.3" case (<fi — 1.89) leads to 
NSE freezeout with appreciably abundant a's compared 
to heavy nuclei (hereafter a-rich QSE; 5^,/lheavy > ijj- 

6 This is equivalent to a "QSE with too few heavy nuclei" com- 
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Fig. 4. — Final abundances as functions of mass number (upper 
panel) for the same tracer particle in Fig. [3] The result of the un- 
changed model is shown in red, and those for the densities divided 
by a factor of 1.3 or multiplied by factors of 2.0, 10 are shown in 
different colors. The position of 48 Ca is indicated by a dashed line. 
The abundance ratios relative to the "p X 10" case are also shown 
in the lower panel. 

In this condition, light nuclei are so abundant that the 
48 Ca abundance continuously decreases. In contrast, for 
the "p x 2" case (<f> = 0.725), the NSE freezes with a 
smaller a abundance that is comparable to that of heavy 
nuclei (hereafter a-poor QSE; rjFhcavy < 10 Such low 
levels of charged light nuclei allow 48 Ca to continuously 
increase during QSE and even after the QSE freezeout 
by Tg ~ 2. In fact, it is the most abundant isotope in 
the end (Fig. ©. The "p x 10" case (</> = 0.145) gives a 
limiting condition of a-poor QSE, resulting in a robust 
48 Ca creation. Figure 0] displays the final abundances 
for the same cases and those relative the "p x 10" case. 
From the lower panel, it is clearly seen that a more a- 
rich QSE tends to shift more abundances from lighter to 
heavier mass numbers. We conclude, therefore, that a- 
po or QSE is the mec hanism for making 48 Ca as suggested 
bv lMever et~aT1 (|1996f ). 

Figure^] (top) shows the result of sensitivity tests with 
various density scaling factors for all the tracer particles. 
Overall, the synthesis of 48 Ca is quit e sensitive to entrop y 
(or <j>) in the range (/> > 1 as shown in lMever et al.l (|1996f) . 
In the explored range, 48 Ca is more robustly produced 

pared to what would be exp ected in NSE with the sam e density 
and temperature described in Meyer et al. (1996, 1998a). 

7 This is equiv alent to a "QSE with to o many heavy nuclei" 
compared to NSE HMever et al.119 96. 1998a), in which QSE cannot 
reduce the number of heavy nuclei and thus 48 Ca survives. 



in the more n-rich environment, Y c < F onuc = 0.417. As 
an example, the particles with Y c w 0.402 in the "p/1.3" 
case have similar 48 Ca abundances to those of the un- 
changed case with Y e ps 0.417, despite larger cj)'s for the 
former (« 1.77) than the latter (« 1.44). This is a con- 
sequenc e of the aan reacti on more rapid in more n-rich 
matter (jMever et al.lll998al) . leading to a faster consump- 
tion of a's. It is important to note, however, that 48 Ca 
is nowhere the most overproduced isotope (i.e., with the 
greatest X/X & ) in the explored range of our unchanged 
model as shown in Figure [5] (bottom). For Y e < 0.417, 
82 Se or 8 Kr are generally the most overproduced iso- 
topes. Only the "p x 2.0" and "p x 10" cases, corre- 
sponding to (j> < 0.7, lead to the production of 48 Ca as 
the most overproduced isotope over a certain range in Y e . 
In summary, the physical conditions of our ECSN model 
lie in the transition region between the q- rich (<j> > 1) 
and a-poor (4> < 1) QSEs (see Fig. 1 in I Meyer et al.l 

DUD. 

4. CONTRIBUTION TO THE GALAXY 

To examine whether ECSNe can be significant produc- 
ers of 48 Ca in the Galaxy, the final abundances are ob- 
tained by mass integration over all the tracer particles. 
Figure[5] (upper- left) gives the isotopic mass fractions (af- 
ter all rad ioactivities hav e decayed) relative to their so- 
lar values (|Loddersll2003l ). i.e., the "production factors", 
as a function of mass number for our unchanged model. 
The "normalization band" between the maximum (352 
for 86 Kr) and one-tenth of that is indicated in yellow with 
the medium marked by a dashed line. Isotopes that fall 
into this band can be accepted as being made in consid- 
erable amounts by ECSNe. A reasonable "flatness" of 
production factors can be seen between A ~ 60 and 90, 
in which many isotopes fall into the normalization band. 
48 Ca (shown by asterisk) is however located slightly be- 
low the normalization band. This implies that ECSNe 
contribute to the production of 48 Ca no more than 10% 
of its solar content. 

As explored in §[3j the nucleosynthesis of 48 Ca is highly 
sensitive to entropy and our unchanged model appears to 
lie in the edge of the a-poor QSE condition. As a sen- 
sitivity test to entropy, the mass-integrated production 
factors for "p x 1.3", "p x 2.0", and "p x 10" are also 
displayed in the upper-right, lower-left, and lower-right 
panels of Figure El respectively. The result of "p x 1.3" 
satisfies the condition of 48 Ca being in the middle of the 
normalization band, which means that the model can 
accounts for one-third of 48 Ca in the solar system. In 
order to fully account for the origin of 48 Ca, one needs 
to increase the densities (that is equivalent to reduce <^'s) 
by a factor of two (Fig. [5J lower- left). The reduced en- 
tropies lead to a-poor QSE conditions in a broad range 
of Y c , in which 48 Ca is copiously produced (Figure El 
top). In fact, it is the most overproduced isotope for 
0-41 < Y c < 0.42 in this case (Figure El bottom). 50 Ti 
and 54 Cr, the isotopes made in a-poor QSE, are also 
enhanced. The change here corresponds to a reduction 
of the entropies from s ~ 14/ce/nuc to ~ llfce/nuc or 
~ 7 fes/nuc in the range Y c < 0.42. This does not seem to 
be an extreme assumption, although the reduced values 
tend to scratch the lower end of the range of entropies 
that are compatible with neutrino-heated ejecta of EC- 
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asterisk in each panel. The normalization band (see the text) is marked in yellow with the medium value indicated by a dashed line in 
each panel. 



SNe or CCSNe. In the limiting case of "p x 10", 48 Ca 
remains the most overproduced isotope, although it is not 
much increased from the "p x 2.0" case. This is a con- 
se quence of robust 48 Ca production for (f> < 1 as proved 
in IMever et al.l (|1996l ). It is interesting to note that the 
res ult in this case r esembles to that of high-density SN la 
by iWoosievI (|1997fl . Isotopes made in a-rich QSE, such 
as 64 Zn, 70 Ge, 74 Se, and 78 > 86 Kr are all depleted and a 
reasonable flatness of production factors cannot be seen. 

The nucleosynthesis of 48 Ca also depends on r oxp . In 
principle, a slower expansion leads to a later freezeout 
from NSE and, consequently, leads to a-poor QSE more 
easily. Such a consequence, however, is not the case in 
our ECSN model. More slowly expanding matter would 
be exposed to neutrino absorptions for a longer time, 
shifting Y e toward higher values as a result of the a-effect. 
Such a Y e increase would compensate the benefits of a 
slower expansion. One may also consider to increase the 
ejecta masses where 48 Ca is richly produced (Y e 0.40- 
0.42). In fact, fine details of the Y e -AM e j distribution 
could well depend on the spatial resolution in the simu- 
lation or on other factors. It is not plausible, however, 
to enhance Ca without others in our unchanged model. 
As found in the lower panel of Figure [21 the most overpro- 
duced isotope in the range Y c < 0.42 is generally 82 Se or 
84 Kr that already lie on the normalization band (Fig. [5j 
upper- left). 



5. IMPLICATIONS 

We studied the nucleosynthesis of 48 Ca using the ther- 
modynamic trajectories of a self-consistent 2D explo- 
sion model of an ECSN (Janka ct al. 2008; Wa naio et al.l 
1201 lh . Appreciable production of 48 Ca was found in 
a-poor QSE conditions. Comparisons of the mass- 
integrated yields with the solar abundances show under- 
production of 48 Ca. This problem would be cured if the 
entropies in the ejecta were somewhat lower than those in 
our original model. We conclude, therefore, that ECSNe 
can be, at least in part, the astrophysical source of 48 Ca. 
The fact that ECSNe can be also associated with the 
origin of many other isotopes (of light tra ns- iron species 
and p otentially of weak r-process species, iWanaio et al.l 



and p 



seems encouraging. 
Our result also has implications for th e interpretation 
of anomalies of 48 Ca found in meteorit es (|Lee et alJll978t 
IMovnier et aHl20Tol : IChen et al]l20Tl . So far, only hy- 
pothetical high-density SNe la have been considered as 
the cause of these anomalies. ECSNe could be, how- 
ever, an additional source of the anomalies of 48 Ca. EC- 
SNe ma y occur at about one-tenth of the rate of norma l 
CCSNe (jlshimaru fc Wanaiolll99flt IWanaio et al.ll2011l) . 
that is, 1-2 events per millennium, which makes them 
much more c ommon than th e rare-class SNe la (once in 
10,000 years. IWoosievI 119971 ). This would be important 
if the composition of the proto-solar system had been 
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affected by a single or a few nearby SNe. 

More work is needed before final conclusions about 
the role of ECSNe as production sites of 48 Ca can be 
drawn. In particular, further improvements of ECSN 
models (e.g., 3D, general relativity, resolution) will be 
important to elucidate whether appreciable amounts of 
n-rich ejecta of ECSNe satisfy the a-poor QSE condition 



needed for the nucleosynthesis of s Ca 
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